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ABSTRACT 

Context. Tidal disruption events are possible sources of temporary nuclear activity in galactic nuclei and can be considered as good 
indicators of the existence of supermassive black holes in the centers of galaxies. 

Aims. A new X-ray source has been detected serendipitously with ROS AT in a PSPC pointed observation of the galaxy cluster A357 1 . 
Given the strong flux decay of the object in subsequent detections, the tidal disruption scenario is investigated as a possible explanation 
of the event. 

Methods. We followed the evolution of the X-ray transient with ROS AT, XMM-Newton and Chandra for a total period of ~13 years. 
We also obtained 7-band optical/NIR photometry with GROND at the ESO/MPI 2.2m telescope. 

Results. We report a very large decay of the X-ray flux of the ROSAT source identified with the galaxy LEDA 095953, a member of 
the cluster Abell 3571. We measured a maximum 0.3-2.4 keV luminosity Log(L x )=42.8 erg s -1 . The high state of the source lasted 
at least 150 ks; afterwards L x declined as~ r 2 . The spectrum of the brightest epoch is consistent with a black body with temperature 
/tr~0.12keV. 

Conclusions. The total energy released by this event in 10 yr was estimated to be A £ >2x 10 50 erg. We interpret this event as a tidal 
disruption of a solar type star by the central supermassive black hole (i.e. ~10 7 M Q ) of the galaxy. 

Key words. X-rays - Galaxies, X-rays: bursts, Galaxies: clusters: individual: A3571, Galaxies: Active, Galaxies: Nuclei - 



1. Introduction 



is growing evidence that most galaxies with a 
host a supermassive black hole (BH) in their cen- 



There 
bulge 

ter (IKormendy & Richston3ll995L iMagorrian et al.lll998l) . The 
most recent theoretical and observational developments in the 
field of galaxy evolution suggest that during their life, galaxies 
experience one or more phases where the BH is active and pow- 
ered by accretion. One of the effects of this phenomenon are the 
so called active galactic nuclei or quasars where accretion occurs 
on time scales longer than 10 6 ~ 7 yr. It is widely accepted that dur- 
ing this period the BH is fed by an accretion disk. On the other 
hand, during the non-active phase of the galaxy life the BHs can 
be briefly powered by star s tidally disrupted d uring close en- 
counters with the BHs (e.g.lFrank & Rees 1 976b. Recent studies 
(see e .g. IWang & MerrittH2004t lEsquei et al.ll2008t iDonlev et all 
2002) pointed out that these events happen with a rate of ~10" 4 - 
10~ 5 galaxy -1 year . These phenomena can generate strong 
flares in the UV-X-ray band as the result of the onset of an ac- 
cretion disk. The X-ray/UV luminosity of the galaxy is therefore 
expected to grow to values close to those of AGNs (i.e. >10 42 erg 
s _1 ). If the star is Sun-like, tidal disruption can only occur if the 
mass of the BHs does not exceed ~10 8 M since, in this case the 
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tidal radius would typically fall within the Schwarzschild radius 
and the star is swallowed as a whole. In the case of rotating BHs, 
the tidal disruption also can occur for highe r BH masses if the 
star a pproaches from a favorable direction (Beloborodo v et ail 
|1992|) . However atmosphere stripping of giant stars is possible 
also for larger BH masses. When assuming a black body radiat- 
ing at the Eddington luminosity, the temperature of the emitting 
region calculated at the last stable orbit (3 R$ , where R$ is the 
Schwarzschild radius) is -30-300 eV for a 10 6 " 7 M G black hole. 
After the disruption, the mass accretion rate, and therefo re the 
bolometric luminosity, may deca y approximately as r 5 ' 3 dReesI 
1990; lEvans & Kochanekl ll989l). even though de viations from 
this law are expected (see e.g. \Lc^ EMm^. 7 candidate 
tidal disruption events have been detected in the X-ray by the 
ROSAT all-sky survey and in the XMM-Newton slew survey 
(see e.glKomossa & Baddl999tlGreiner et al.l2000llGrupe et al] 
ll999UKomossaetalJl2004HEsquei et al.ll2007l) . Candidate tidal 
disruption have been observed also in the Optical/UV (e.g. 
iRenzini et alJll995t iGezari et al.l l2008). The X-ray events show 
strong X-ray variability (i.e. more than a factor of 100), 0.3- 
2.4 keV luminosities typical of AGN and quasars (i.e. >10 42 erg 
s ) declining as ~r 5 ^ 3 , ultrasoft (i.e. ~0.1 keV) X-ray spectra 
and the absence of Seyfert activity in ground-based optical spec- 
tra (see Komossa 2002 for a review). In this letter we report the 
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Fig. 1. Left Panel: The ROSAT-PSPC 0.3-2.4 keV colour coded 
image of the field of A3571. The X-ray transient is indicated 
with the arrow. Right Panel: GROND image in the K,$ -band of 
the the region of the X-ray transient, the red circle represents 
the ROSAT-PSPC error-box and the red ellipsoid is the Chandra 
confidence region. 
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Fig. 2. Left Panel: The ROSAT-PSPC unfolded spectrum of 
the source TDXFJ 134730.3-325451 (Crosses) and the best fit 
blackbody model Continuous Line. 



serendipitous detection of a very large d ecay in the X-ray flux 
of the galaxy LEDA 095953 at z=0.0366 (lOuintana & de Souzal 
[1991) in the field of view of the galaxy cluster Abell 3571 us- 
ing ROSAT, XMM-Newton and Chandra archival X-ray data and 
GROND optical/NIR follow-up. 

2. Observations and data processing 

During the soft X-ray follow-up of t he Swift-BAT, har d X- 
ray selected sample of galaxy clusters dAiello et alj|2008l) . we 
serendipitously discovered a bright source in a ROSAT-PSPC 
image of Abell 3571 (see Fig.Q]) that apparently disappeared in 
subsequent observations with XMM-Newton and Chandra. The 
field of Abell 357 1 has been observed in X-rays in pointed mode 
with the ROSAT-PSPC/HK1 detectors, XMM-Newton EPIC, 
Chandra ACIS-S and during the ROSAT all-sky survey. The 
ROSAT data were processed using the standard MIDAS-EXSAS 
software version 03SEPpll, while XMM-Newton and Chandra 
data were processed using XMM-SAS version 7.0 and CIAO 
4.0. We also performed an optical/NIR follow-up of the field 
in the g', r, i', z' ba nds and J, H, Ks bands, by using GROND 
dGreiner et al.ll2008b . mounted at the 2.2 m ESO/MPI telescope 
at La Silla observatory (Chile). In Table Q] we list the log of the 
observations of the field of Abell 3571 relevant for this work. 



2.1. Data analysis 

ROSAT-PSPC: Using the ROSAT-PSCP 0.3-2.4 keV energy 
band image and exposure map available in the archive, we cre- 
ated a background map, by first running a sliding cell source 
detection with a signal-to-noise ratio threshold of 4, excising 
the detected sources from the image and fittingthe residual data 



Instrument 


Date 




Band 


Exposure 




dd-mm-yyyy 






s 


ROSAT-PSPC 


12-08-1992 


0.3 


-2.4 keV 


5800 


ROSAT-HPJ 


05-08-1994 


0.3 


-2.4 keV 


19200 


XMM-Newton-EPlC-PN 


29-07-2002 


0.3 


-2.4 keV 


24200 


Chandm-AClS-S 


31-07-2003 


0.3 


-2.4 keV 


33900 


ESO/MPI-GROND 


12-08-2008 


g. 


r',i\ z' 


458 


ESO/MPI-GROND 


12-08-2008 




J,H,K 


480 



Table 1. The log of the X-ray and optical/NIR observation of 
the field of Abell 3571 used for this work. 



with a spline model. With the background map described above 
and using the exposure map we run a maximun likelihood, PSF- 
fitting, source detection. The bright, eye-detected, source men- 
tioned in the previous section has been detected with a signif- 
icance >10cr with coordinates of the centroid 13 /l 47 m 29.8 s - 
32°55'00" (J2000) with an error box of -25". Considering that 
the source is embedded in the X-ray emission of the galaxy 
cluster Abell 3571, the photometric properties have been eval- 
uated using aperture photometry instead of using the estimate of 
the detection software. We extracted the source counts in a cir- 
cular region with a radius corresponding to the 90% Encircled 
Energy Radius (EER) around the centroid estimated by the max- 
imum likelyhood fit. Background counts have been estimated in 
a ring around the source with an inner radius 1.5 times the 90% 
EER and outer radius 2 times the 90% ERR. As a result we ob- 
tained 305+27 net source counts corresponding to a count-rate 
of (5.2+0.4)xl0~ 2 s -1 . PSF fitting estimated a flux -30% lower 
as a result of a likely overestimate of the background. In the same 
regions we also extracted source and background lightcurves and 
spectra. All the epochs but one have flux consistent with a mean 
count-rate of ~5.2xl0~ 2 s _1 . Note that the observation was split 
in two parts separated by ~150 ks. This indicates that no strong 
variability is present on such a time scale. However in one time 
bin the count-rate reached (8.3+1.0)xl0~ 2 s _1 . This flux is ~3<x 
above the mean and could indicate a small flare during the event. 
The unfolded spectrum of the source is shown in Fig. [2] . 
We performed a x 2 fit to the spectrum using a black-body 
model plus cold absorption. We kept as free parameters the 
temperature kT, the columns density N# and the normaliza- 
tion k. As a result we obtained a temperature kT=120+16 eV 
and N H =3.93+! ^xl0 20 cnr 2 and * 2 /v= 1.05/4.00. This corre- 
sponds to a 0.3-2.4 keV flux of 7.65+0.92x10"° erg ctrT 2 s _1 . 
Assuming no spectral changes during the observation, the maxi- 
mum 0.3-2.4 keV flux emitted by the source is 1.16±0.13xl0 -12 
erg ctrr 2 s _1 . Note that the column density is consistent within 
lcr w ith the Galactic value predicted by Dic key & Lock man 
( 1990) (i.e. N#=3.7xl0 20 ctrT 2 ) suggesting an extragalactic ori- 
gin of the source. We also tested other spectral models, like a 
power-law plus absorption, without improving the quality of the 
fit by having F ~6.5, and N# ~3xl0 21 cm" 2 as best fit parame- 
ters and * 2 /v=5. 30/4.00. 

We also performed aperture photometry on ROSAT all sky sur- 
vey images at the source location and we determined a lcr upper- 
limit (see next Section) on the count-rate of 0.44 s _1 correspond- 
ing to a flux of ~2.2xl0~ 12 erg cm - 
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assuming as a spectral model a black body with kT=0.12 keV. 
ROSAT-HRI: The source is visually observable also in the 
HRI observation though much fainter than in the PSPC image. 
Considering the highest spatial resolution of the HRI detector 
we performed a source detection to constrain the source position 
better than in the PSPC pointing. Unfortunately because of the 
faintness of the source and the difficulty in modeling the high 
background due to the extended X-ray emission of the cluster, 
the source has been detected with a significance of ~3-4<x, but 
as an extended source. We therefore performed aperture photom- 
etry around the ROSAT-PSPC centroid. By adopting the same 
procedure used for the PSPC data, we measured 38+9 source 
counts corresponding to a count rate of (1.97+0.45)xl0~ 3 s _1 . 
By assuming no spectral change between the ROSAT-PSPC and 
HRI observation (i.e. 2 years) we estimated a 0.3-2.4 keV flux 
of 1 .53+0.4X 1(T 13 erg ctrr 2 s" 1 . 

XMM-Newton: A visual inspection of the EPIC-PN image 
does not show any relevant signature of a X-ray source in the 
ROSAT error box. We run a maximum likelihood source detec- 
tion in that region but no sources were detected with signifi- 
cance >3cr. We therefore estimated the 3<x upper-limit of the 
so urce count at t he source location by using the prescriptions 
of lNarskvl (l2000h. Given M coun ts actually observed in a re- 
gion of 30" ICappelluti etai1l2009l) . and B background counts, 
the lcr upper limit is defined as the number of counts X that 
gives the probability to observe M (or fewer) counts equal 
to the formal 68.3% Gaussian probability: P(< M,X + B) = 
PcausA^- 3%). Assuming Poissonian statistics this equation be- 
comes: P Causs = <r (z+B) Yjfio ^tt^- By iteratively solving the 
previous equation in the case of Pc m <ss=0-997, we obtained the 
3<r upper limit X=167 counts. This corresponds to a count-rate 
of 8.8xl(T 3 s _1 and to a flux of 1.45xl0" 14 erg ctrT 2 s^Q 
Chandra: Using the task wavdetect, included in the CIAO soft- 
ware, we performed a source detection on the ACIS-S4 Chip 
where the source is observable. The analysis was performed 
in the 0.3-2.4 keV energy band. We ran wavdetect using an 
exposure map estimated at 1 keV with a threshold of 10~ 6 , 
corresponding to a maximum of 1 spurious detection expected 
per ACIS-chip. As a result, in the ROSAT PSPC error box we 
detected one source with 9.5+3.74 counts corresponding to a 
count-rate of 2.89xl0~ 4 s _1 . The new source coordinates are 
tt=13*47 m 30.33*, <5=-32°54'50.63". At the source position the 
estimated cr positional error is ~ 2.5". We also ran a source 
detection in the 2-7 keV band without detecting a source. We 
also estimated the flux of the source by folding into XSPEC 
the Chandra response matrices and assuming as the spectral 
model the ROSAT-PSPC best fit and obtained a 0.3-2.4 keV 
flux of 3.43±1.35xl0~ 15 erg crrT 2 s -1 . Note that the source cen- 
troid lies about 5" from the chip edge, therefore the estimate 
of the background could be affected by features near this edge. 
According to its coordinate, the X-ray source has been named 
TDXFJ 134730.3-325451. 

GROND: In Fig. □ we show the GROND K-band image 
of the area covered by the ROSAT-PSPC error box. As one 
can notice, the Chandra ellipsoid is placed right on the disk 
of a bright galaxy with an offset of ~3 " from the op- 
tical centroid of the galaxy. Note that the GROND images 
have an astrometric uncertainty of 0.5". We searched the 
NED catalog and the galaxy was previously known as LEDA 
095953, a member of the cluster Abell 3571 at z=0.0366 




10" 

waveleng th A [A] 

Fig. 3. The SED of the galaxy LEDA 095953 in the g', r', i', z', 
H, J and K s bands obtained with GROND over the best fit SED 
model of a SO galaxy {solid line). 



dOuintana & de Souzal Il993l) . We determined the photomet- 
ric properties of the source and obtained the following AB 
magnitudes: 77.05+0.10, 16.26+0.10, 15.80+0.10, 15.54+0.10, 
15.04+0.02, 14.75+0.01 and 14.87+0.04 in the g' ,r', i', z', 
J, H, bands, respectively. These magnitudes are calculated 
by using SA105-815 and 2MASS stars as a reference and cor- 
rected for the expec ted foreground extinction of £Vy=0.054 
dSchlegel et aljfl998l) . We examined the global spectral energy 
distribution (SED) in Fig.[3]of this object using a tool based on 
a chi-squared minimization method that allows us to fit the ob- 
served fluxes with a combination of AGN and galaxy emission, 
also allowing for the possibility of extinction of the AGN flux. 
From this analysis we found that the model fit that better de- 
scribes the SED of this object is a normal SO galaxy with no 
significant contribution due to the presence of an AGN. 



3. Long term light curve 

The X-ray flux of TDXFJ134730.3-325451 has been followed 
for 13 years and its lightcurve is presented in Fig. [4] Remarkably, 
it shows an amplitude of the variability of a factor >650. We esti- 
mated the maximum 0.3-2.4 keV luminosit>0 as Log(L z )=42.83 
erg s _1 . By excluding the upper limits, we fitted the long term 
lightcurve with a power-law model in the form: Lx(t) — k * 
('-jff-y a , where k is the normalization, t<n sr marks the time of 
the star's disruption and a is the slope of the decay. As a re- 
sult we obtained k~7xl0 42 erg s" 1 , ^, sr ~1991.6 and a=2. 2+0.5. 
Because of the low number of data points, the uncertainties on 
the fit parameters are large. The value of the slope a is therefore 
marginally consistent with the predicti ons of variants of tidal dis- 
ruption m odels (i.e. a ~5/3, see e.g. lEvans & Kochanekl 1989: 
iRee si [19901) . Note that according to the lightcurve, the source 
at the maximum could have been at 1-2 orders of magnitude 
brighter with respect to the time of the ROSAT-PSPC observa- 
tion. If the source continues its decay, it would be detectable 
with Chandra for all the year 2009 and 2010, with moderately 
deep exposures. 



1 The flux was derived from the count-rate with the PIMMS 
software by assuming as the spectrum the .ROSAT-PSPC best fit. 
http://heasarc.nasa.gov/Tools/w3pimms.html 



2 The luminosity has been computed assuming a A-dominated uni- 
verse with n A =0.7, n„,=0.3 and H =72 km s~' Mpc _1 
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Fig.4. The X-ray lightcurve of TDXFJ134730.3-325451. The 
solid curve represents the result of the power-law fit while the 
vertical line marks the expected begin of the flare. 



4. Discussion 

We serendipitously detected a factor >650 decay in the flux of 
an X-ray source in the rich galaxy cluster A357 1 . The flare most 
likely originated in the galaxy LEDA 095953 at z=0.0366. The 
most luminous epoch in the lightcurve reached LogLx=42.83 
(0.3-2.4) and lasted for more than 150 ks, although higher lumi- 
nosities at the time of the burst cannot be excluded. The spec- 
trum of the flare is consistent with a black-body of temperature 
kT~0.12 keV with no significant hard components. The opti- 
cal/NIR SED of the host galaxy does not show any signature 
of AGN activity. In addition to this, the decay of the long term 
light curve is marginally consistent with t~ 5 ^ 3 , and thanks to 
the stringent RASS upper limit, we likely detected the source 
close to the burst. The X-ray luminosity and the variability of 
the source are therefore much higher than those of ULX, which 
have Lx <10 39 " 40 erg s" 1 and flux variation s typically of up to a 
factor of 2-3 (see e.g. Colbert & Ptak 2002). Supernovae explo- 
sions can also produce X-ray flares, but these are either of very 
short duration (hours), or else if longer lasting they have much 



harder X-ray spectra an d typically have Lx<10 41 erg s _1 (see e.g. 
Imml er & Lewinl2003l for a review). Additionally novae and su- 
pers oft sources can produce bur sts but with much shorter dura- 
tion dTanaka & Shibazakilll996l) . According to the event prop- 
erties, we can safely exclude these kinds of stellar phenomena. 
However, such a strong variability could be observed in GRB 
afterglows or Hypernovae but with a much shorter duration of 
the outburst. The characteristics of the flare are therefore consis- 
tent wUhj)reviOTi^cis^ events (see 
e.g. iKomossa et al" 20041: lEsquei et al.1l2008l) a nd with theoreti- 
cal predictions dReeslll990HLodato et alj|2008l) . The mass of the 
BH should be not much higher than ~10 8 M to disrupt a Sun- 
like star outside the Schwarzschild radius and to allow the onset 
of accretion. We then estimated the black hole mass fr om the 
Lk-Mbh relation determined by Marc oni & Hunt! d2003l) in the 
local Universe. If is the K-band luminosity of the bulge of the 
galaxy in solar units then log(M BH /M B ) -8.08+ 1.2 1( L Kbiilge [L Q - 
10.9). Using an aperture radius of ~4" (i.e. ~2.5 kpc, at the 
galaxy distance) we estimated a log(L/s- M!f )~10L/L o which cor- 
responds to log( M bh I 'M Q )~1 . The total energy released during 

r-CO 

the event can be estimated with AEx-j l Lx(t)dt, where Lx(t) 
is the best fit to the light curve. By integrating over a period of 
~10 yr after our first data point we obtained a total energy re- 



lease of the order of 2xl0 50 erg. If e~0.1 is the typical efficiency 
of mass to energy conversion for accretion onto a BH, then the 
total mass accreted during a 10 year period is AM~ AEx/ec 2 . 
This yields a total mass deposition of lxlO" 3 M . This value is 
actually a lower limit, since we started the integration only at the 
start of the observations (note that the fitting function diverges 
at the time of the disruption). Moreover it is possible that the 
total emission is dominated by an unobserved EUV component 
(extrapolation of the single black body fit itself predicts a 7% 
higher luminosity), and the t rue accretion rate is much higher. 
According to the results of iLi et al (2002), for a black hole of 
10 7 M Q , the flare occurs ~0.07 yr after the disruption. A 10 yr 
integration of the lightcurve model from this epoch yields an 
estimate o f the accreted ma ss of ~0.03 M . A numerical sim- 
ulation by lAval et al] ((2000) showed that for a 1 M star dis- 
rupted by a 10 6 M Q BH, only a fraction of the order of ~1Q% of 
the initial mass is actually accreted. Our limits on the accreated 
mass are therefore consistent with the scenario of a disruption 
of a 1 M star by the BH. However a partial stellar disruption 
or an explosive disruption dBrassart & L uminet 2008) cannot be 
ruled out, although disruption and accretion of a brown dwarf 
or of a planets is very unlikely. We note that since the tidal 
disruption rates are likely of the order of 10~ 4 -10" 5 galaxy" 1 
yr" 1 , galaxy clusters are ideal laboratories to explore this phe- 
nomen on: by conservatively assumin g 100 member galaxies per 
cluster dReiprich & Bohringeril2002l) . one should expect 1 tidal 
disruption e vent per cluster ever y 10-100 yrs. The eROSITA all- 
sky survey dPredehl et al.ll2006h will detect several thousands of 
galaxy clusters and therefore a large number of these events. 
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